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Abstract 

We study the effect of the spreading widths on the excitation prob- 
abilities of the double giant dipole resonance. We solve the coupled- 
channels equations for the excitation of the giant dipole resonance 
and the double giant dipole resonance. Taking Pb + Pb collisions as 
example, we study the resulting effect on the excitation amplitudes, 
and cross sections as a function of the width of the states and of the 
bombarding energy. 

Double giant dipole resonances have been mainly studied in heavy ion 
Coulomb excitation experiments at high energies (for a recent review, see 
0). The feasibility of such experiments has been predicted in 1986 || || 
where the magnitude of the cross sections for the excitation of the Double 
Giant Dipole Resonance (DGDR) was calculated (see also In ref. || a 
recipe was given for treating the effect of the width of the giant resonances 
on the excitation probabilities and cross sections. In this letter we make a 
quantitative prediction of this effect using a realistic coupled-channels calcu- 
lation for the excitation amplitudes. The coupling interaction for the nuclear 
excitation i — > / in a semiclassical calculation for a electric (n = E), or 
magnetic (ir = M), multipolarity, is given by (eqs. (6-7) of ref. |J) 
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n-,v (» = i I ) ' ^ - ^J^p^Q^ (^r)M(nX^) . (2) 

Above, 6 is the impact parameter, 7 = (1 — /5 2 ) , /3 = v/c, r = <yvt/b 
is a dimensionless time variable, eo = '"fhv/b sets the energy scale and 
Qir\/j.(£, T ), with £ = &f = (Ef — Ei) /eo as an adiabatic parameter, de- 
pends exclusively on the properties of the projectile-target relative motion. 
The multipole operators, which act on the intrinsic degrees of freedom are, 
as usual, 

M(E\,fi) = J d 3 r p(r) r x Y^r) , (3) 

and 

M(Ml,/i) = ~ J dh J(r).L(rF lM ) , (4) 

We treat the excitation problem by the method of Alder and Winther||. 
We solve a time- dependent Schrodinger equation for the intrinsic degrees 
of freedom in which the time dependence arises from the projectile-target 
motion, approximated by the classical trajectory. For relativistic energies, 
a straight line trajectory is a good approximation. We expand the wave 
function in the set {| k); k = 0, N} of eigenstates of the nuclear Hamiltonian, 
where denotes the ground state and N is the number of intrinsic excited 
states included in the coupled-channels (CC) problem. We obtain a set of 
coupled equations. 

To simplify the expression we introduce the dimensionless parameter 0jy 
by the relation 
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Then we write the coupled channels equations in the form 

N 



da k (r) 
dr 



T) CLAT) . 



(6) 



j=0 irXfi 



In what follows we concentrate on the El excitation mode. In this case, 
we have 



r^(r)-^(^) 2 0(r) 



Qei±i{^t) = T ^{t) 



(7) 

where </>(r) = (1 + r 2 ) 1/2 . 

Following ref. || the inclusion of damping leads to the coupled-channels 
equations with damping terms, i.e., 



da k (r) 
dr 



N 



^QEi^kj,?)®^ exp(z^r) aj{r) -a k (r). (8) 
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These equations lead to master equations for the occupation probabilities, 
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Pj(r) = \aj(r)\ , in the form 

dP k {r)/dT = G k {r)-L k {r) 

where 



(9) 



G k (r) = 2 3m ^ ^ Q^kj, r)e^f exp(^r) a k (r) a*(r) (10) 
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and 



L k (r) = -—P k (r) . 
These equations can be integrated, yielding the conservation law, 



(11) 
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where Fu(r) 



^ / P k (r')dr' (12) 



Due to the exponential decay of the states with k ^ 1, we have for t — > oo 
the limit -Pfc(oo) = <5j P (oo) and 

P (oo) + ^P fc (oo) = l. (13) 

This means that for i — > oo there is a probability to find the system in 
the ground state given by Po(oo) an d a probability that it has been excited 
and decayed through the channel j which is given by Fj (oo). Thus, the set of 
equations 13 are shown to correctly describe the contribution to the excitation 
through channel j. 

The excitation probability of an intrinsic state | j) in a collision with 
impact parameter b is obtained from an average over the initial orientation 
and a sum over the final orientation of the nucleus, as 

W) = 9fVr £ l*k°°)l a > ( 14 ) 

and the cross section is obtained by the classical expression 

Oj = 2tt j P^b) T(b) bdb . (15) 

Above, Tib) accounts for absorption according to the prescription of ref. ||, 
using the nucleon-nucleon cross sections and the ground state density of Pb 
from experimental data. 

We consider the excitation of giant resonances in 208 Pb projectiles, inci- 
dent on 208 Pb targets at 640 A-MeV. This reaction has been studied at the 
GSI/SIS, Darmstadt @. For this system the excitation probabilities of the 
isovector giant dipole (GDR) at 13.5 MeV are large and, consequently, high 
order effects of channel coupling should be relevant. To assess the importance 
of the damping effects, we calculate the matrix elements assuming that the 
GDR is an isolated state depleting 100% of the energy-weighted sum-rule. 
The matrix element for the GDR— > DGDR transition incorporates the boson 
factor v^2, as usual j3]. The energy location of the DGDR state is taken as 
27 MeV, consistent with the experimental data. The spin and parities of the 
states are given by 1~ for the GDR, and + and 2 + for the DGDR, respec- 
tively. The distribution of the strength among the + and 2 + DGDR states 
are simply obtained from Clebsh-Gordan coefficients M. 
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In figure 1 we plot the time-dependent occupation probabilities of the 
ground state, in Pb, N = 0, of the GDR state, N = 1, and of the DGDR 
state, N = 2, respectively. Figure 1(a) shows the occupation probabilities 
with the widths equal to zero, = 0. In figure 1(b) we plot the occupation 
probabilities of the GDR state, N — 1, and of the DGDR state, N = 2, 
with T N = (full lines), and with T GDR = 4 MeV (experimental), and 
Tdgdr = 5.7 MeV (dashed lines). The width of the DGDR is set to T DGDR = 
V2Tcdr, following the apparent trend of the experimental data 0. Note, 
that Tdgdr = ^gdr has a better (and simpler) theoretical explanation ||. 
But, the time integrated population of the DGDR state will not be much 
influenced by using the later parametrization. 

We observe that the inclusion of damping leads to strong modifications in 
the time-dependent occupation probabilities of the GDR and DGDR states. 
One might wrongly deduce from figure 1 that the excitation probabilities 
are reduced proportionally to the difference between the maximum value of 
Pn(t), with and without damping. However, the quantities shown in figure 
1(b) includes the loss of occupation probabilities of a given state, while it 
is being populated by the time-dependent transitions. Thus, the reduction 
of the excitation probabilities of the GDR and the DGDR due to damping 
is smaller than deduced from figure 1. The relevant quantity to calculate 
the excitation cross sections are the quantities Fj(oo), which account for the 
time-integrated transition probability, followed by decay, of the state j. 

In figure 2 we plot the flux functions, or time-integrated transition prob- 
abilities to the GDR and the DGDR states as a function of the width of 
the collective state, Tgdr, keeping constant the ratio Tdgdr/^gdr, = \/2. 
We keep the impact parameter fixed, b = 15 fm. We note that, varying the 
width from up to 4 MeV leads to a 10% decrease of the flux functions 
into GDR and DGDR states. A similar tendency is observed for the total 
cross section, integrated over impact parameters. The excitation probability 
decreases with about l/b 4 , therefore the grazing collisions are weighted most 
strongly. 

In figure 3(a) we plot the effect of damping in the total cross sections, as 
a function of the bombarding energy. We take Tgdr = (dashed line) and 
Tgdr = 4 MeV (full line), keeping constant the ratio T DGDR /T GDR = \/2. 
We observe that the effect of damping disappears, as the bombarding energy 
increases. At high energies the reaction is fast, and the system does not have 
time to dissipate during its excitation. In this regime the sudden approxi- 
mation is a valid approach to the calculation of the excitation amplitudes. 
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We note that our model is restricted to the excitation of isolated res- 
onant states, including a time- dependent loss term on the far right of the 
coupled-channels equations |8|. This is different from a study of the influence 
of the fragmentation of the resonances into many neighbouring states. In 
this case, the Coulomb excitation of the giant resonances is obtained as a 
superposition of excitations to states spread over an energy envelope, usually 
taken as a Lorentzian shape. States at lower energy are more easily excited 
than states at higher energies. Thus, the spreading of the resonances may 
lead to another kind of effect of the widths, not obtainable in the above 
treatment. To study this effect in a simple way, we use the harmonic model 



of ref. |T(|. The photo- nuclear cross sections which enter in these calcu- 
lations are of Lorentzian shapes, with a cut at the low energy limit of 8 
MeV, corresponding to the threshold for neutron emission, since most of 
the DGDR manifestation in experiments come from neutron emission after 
relativistic Coulomb excitation. The magnitude of the photo-nuclear cross 
sections are obtained by using a 100% depletion of the Thomas-Reiche-Kuhn 
energy- weighted sum-rule applied to the GDR. 

The harmonic model provides a simple analytical formula to calculate the 
excitation probabilities of the DGDR [|K|. The resulting cross sections are 
shown in figure 3(b) where in the solid line we take Tqdr = ^dgdr = 0, while 
the dashed-lines are for T DGDR = 5.7 MeV and T GDR = 4 MeV. We observe 
a similar effect as in figure 3(a). This is understood as a reduction due to 
the spread of states at energies above the energy centroid of the Lorentzian 
envelope. The excitation amplitude for these states are smaller, thus leading 
to a net reduction of the energy integrated Coulomb excitation cross sections. 

In conclusion, we have obtained the dependence of the excitation ampli- 
tudes on the width of the giant resonance states. We show that the effect 
reduces excitation probabilities, and cross sections. We have developed an 
approach to solve this problem in realistic situations. It is demonstrated that 
the dynamical effect of the widths of the GR's in a time-dependent picture 
leads to a decrease of the cross sections, more accentuated for low energy col- 
lisions. The energy fragmentation of the giant resonances can be studied in a 
simple fashion within the harmonic model. The net effect is also to decrease 
the cross sections with increasing width, specially at low energy collisions. 
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Figure Captions 

Fig. 1 - Occupation probabilities in Coulomb excitation of Pb, in Pb + Pb 
collisions at 640 MeV.A. N = for the ground state, AT = 1 for the GDR 
state, N = 2, for the DGDR state, respectively. Figure 1(a) shows the 
occupation probabilities with the widths equal to zero, = 0. In figure 
1(b) we plot the occupation probabilities of the GDR state, N — 1, and of 
the DGDR state, N = 2, with T N = (full lines), and with T GDR = 4 MeV, 
and T DGDR = 5.7 MeV (dashed lines). The width of the DGDR is set to 
Fdgdr = Vz^cdr, according to the trend of the experimental data [0]. 

Fig. 2 - Flux functions, or time-integrated transition probabilities to the 
GDR and the DGDR states in Pb as a function of the width of the collective 
state, Tgdr, keeping constant the ratio r ' dgdr/^gdr = v2- We keep the 
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impact parameter fixed, b = 15 fm. N = 1 for the GDR state, N = 2, for 
the DGDR state. 

Fig. 3 - (a) Total cross sections, as a function of the bombarding energy. 
We take Tgdr = (dashed line) and T GR)R = 4 MeV (full line), keeping 
constant the ratio Tdgdr/^gdr = V%- N = 1 for the GDR state, N = 2, 
for the DGDR state, (b) Results of the harmonic model [10 for the Coulomb 
excitation cross sections of the GDR (N = 1) and the DGDR (N = 2). The 
solid lines use T GDR = T DGDR = 0, while the dashed-lines are for T DGDR = 
5.7 MeV and T GDR = 4 MeV. 
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